Dimensions, Maximal Growth Sites and Optimization in the Dielectric Breakdown 

Model 



Joachim Mathiesen 1 , Mogens H. Jensen 2 , and Jan 0ystein Haavig Bakke 3 
1 Physics of Geological Processes, University of Oslo, Oslo, Norway 
2 The Niels Bohr Institute, Blegdamsvej, Copenhagen, Denmark 
3 Department of physics, Norwegian University of Science and Technology, Trondheim, Norwa^\ 

We study the growth of fractal clusters in the Dielectric Breakdown Model (DBM) by means of 
iterated conformal mappings. In particular we investigate the fractal dimension and the maximal 
growth site (measured by the Hoelder exponent a m in) as a function of the growth exponent n of the 
DBM model. We do not find evidence for a phase transition from fractal to non-fractal growth for a 
finite Tj-value. Simultaneously, we observe that the limit of non-fractal growth (D — > 1) is consistent 
with Oimin — * 1/2. Finally, using an optimization principle, we give a recipe on how to estimate the 
effective value of n from temporal growth data of fractal aggregates. 



I. INTRODUCTION 

Laplacian growth and the formation of complex patterns 
has been the subject of numerous theoretical and exper- 
imental works. The classical examples are the ramified 
pattern appearing in a Hele-Shaw cell when a less vis- 
cous fluid is injected into a more viscous fluid |l| and 
the fractal structures emerging from the particle aggre- 
gation in Diffusion- limited Aggregation (DLA) Q. In 
the latter example mono-disperse particles are released 
one-by-one from a remote source and diffuse until they 
hit and irreversibly adhere to a seed cluster at the center 
of coordinates. The cluster slowly expands as particles 
are added. Statistically, the motion of a single particle 
is described by the harmonic potential U satisfying the 
Laplace equation AU = and the probability of sticking 
to the cluster at a specific site, z, is given by the harmonic 
measure |Vt/(z)|. The formulation of DLA is contained 
within a more general model, the Dielectric Breakdown 
Model (DBM) [g, 0] , where the growth probability p v at 
the cluster interface is proportional to the harmonic mea- 
sure raised to a power 77, p v cx |V£/|''. Despite intensive 
research in Laplacian growth, fundamental questions re- 
garding the scaling properties still have no answer. The 
growth laws of DLA and DBM are extremely simple and 
in apparent disparity to the complex patterns they pro- 
duce. The complex patterns arise from a strong correla- 
tion between the position of already aggregated particles 
and the influx of new particles. As the outermost tips ad- 
vance the probability for particles to reach the parts left 
behind diminishes and the harmonic measure broadens 
and becomes even multifractal [5|. For increasing values 
of T) the growth probability will concentrate around the 
tips and the fractal dimension gets closer to unity and 
ultimately, in the limit of infinite 77, the particle cluster 
loses fractality. Recently, it has been speculated that in 
two dimensions this transition from fractal to non-fractal 
growth may happen at a finite critical value of 77 and nu- 



merically, this value has been found to be r\ fts 4 0, @] . 
In the vicinity of such a critical point it may be safe to 
disregard the noise giving rise to local density fluctua- 
tions along the branches 0]. For that reason, the domi- 
nating stochastic component in the cluster growth is the 
rate at which growing tips split in two or more branches. 
While growing, neighboring branches compete and if one 
branch quickly dies after a tip-splitting the growth will 
stay non-fractal. It has been shown [8| that in the ide- 
alized case of straight growing branches, tip-splitting is 
suppressed for 77 > 4 supporting that r\ c = 4. Based 
on the idealized branch growth model a renormalization 
group approach has been used in an expansion around r; c . 
Although an expansion provides important information 
for small values of 4 — 77 it may provide little information 
on DLA (77 = 1). In this article we test the hypothesis of 
a critical point at 77 = 4 performing extensive numerical 
simulations. We provide detailed figures on the depen- 
dence of the fractal dimension, a m i n and the exponent 
r\. Moreover, we propose a method for extracting effec- 
tive 77-exponents given either experimental or numerical 
data series. For that purpose we make use of iterated 
conformal maps Q which have proven a convenient tool 
for generating conformal mappings of domains of arbi- 
trary shape [l0( , see section [TTJ In Section IIIII a method 
is proposed for extracting effective r\ exponents by opti- 
mization. In section IIVI we present results pro et con a 
phase transition in DBM, the maximal growth sites and 
the fractal dimensions. 



II. ITERATED CONFORMAL MAPPINGS 

The conformal invariance of the Laplace equation re- 
duces the problem of finding the harmonic measure, 
pi(z), around any simply connected domain in the com- 
plex plane to that of finding a conformal transformation 
u> = < J > ~ 1 (z) of the domain to the unit disc 
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The method of iterated conformal mappings provides a 
general framework to construct such transforms as well as 
a simple procedure to grow DLA clusters. Assume that 
a DLA cluster of n particles is mapped to the unit disc 
by < &~ 1 . An extra particle is added to the cluster by first 
adding a small bump of size \/X to the unit disc using 
a mapping y n +i and subsequently applying the inverse 
mapping Finally, the composed mapping $„ o ip n+ i 
transforms the unit disc into a cluster of n + 1 particles. 
The basic mapping tp n +\ is defined by two parameters 
the position and size of the bump, the position, e 10 , is 
random in DLA since the measure is uniform around the 
circle. The size y/\ of the n'th bump is controlled by the 
condition that 

Va^ = \/a„|$v ;9 )| 

Consequently, the particles (transformed bumps) will all 
to linear order have the same size V^o- The full recursive 
dynamics is written as iterations of the basic map 



<5> (n) (V) = (fe u Xi ° ■ ■ ■ ° fe n ,\ n (w) 



(2) 



Note that this structure is unusual in the sense that the 
order of iterates is inverted compared to standard dy- 
namical systems. 

For DBM the growth measure along the cluster interface, 
parameterized by s, is given by 



Pv( s ) 



(3) 



which for n ^= 1 is not conformally invariant. On the unit 
circle, parameterized by 9, the growth measure trans- 
forms into 



p v (6)dd ~ p v (s(6)) 



ds 
Iff 



de~ |$V )r~ ? 'd6> . (4) 



In the simulations we choose 9 according to the distri- 
bution p n using standard Monte Carlo samplings of the 
measure. The number of samples needed for an accu- 
rate estimate of the distribution increases with n and is 
chosen according to 



(5) 



/A^max s pi(s) 



By choosing k > 1, the site of maximal measure will on 
the average be visited more than once during the sam- 
pling. It turns out that there is no visible change in the 
scaling of the clusters when choosing k > 1, see Fig. Q] 
for a test of convergence as function of fc; in the results 
presented below, we use 2 < k < 8. 



III. EXTRACTING EFFECTIVE rj EXPONENTS 
BY OPTIMIZATION 



Consider an interface growing at a rate determined by 
some unknown function of the harmonic measure. The 
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FIG. 1: (Color online) Test of convergence for the deployed 
Monte Carlo method. Estimates of the fractal dimension for 
rj — 4 and for Monte Carlo samples given by k — 1/4, 2, 8, 32 
in Eq. Q. For each data point, we used four clusters of size 
40000 particles. 



method of iterated conformal mappings is readily turned 
into a framework for estimating this function. More 
specifically, it is here demonstrated on numerical simula- 
tion data of the DBM that the value of n can be extracted 
from a careful tracking of the cluster growth. The gen- 
eral idea is to utilize the iteration scheme in tracking the 
motion of the interface by gradually expanding the map- 
ping, see [Io| for further details. The harmonic measure 
is recorded as the interface evolves and from a maximum 
likelihood principle the n value of the growth is extracted. 
The probability for growth to occur at a site z n at the 
interface is in a given growth step n approximated by the 
sum 



Pn(n,z n ) 



p\(n,z n ) 



(G) 



From this expression, more ways exist to estimate the 
rj value used in the simulation. Assuming that the n'th 
growth event occurred at the site z n , a direct estimate 
of rj follows from maximizing p v (z n ) with respect to rj. 
Naturally, this will lead to dramatic fluctuations in the 
estimates and therefore maximizing products of p^ over 
several growth steps provides a better estimate, 



Y[Pr,(k,Z k ) 



(7) 



In Fig. [21 we show how this product varies as function 
of rj and with the number of factors used. With an in- 
creasing number of factors the maximum becomes more 
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FIG. 2: Product of the growth measure, maximum likelihood, 
U.k=i Pr,(?k,z 2 k) as function of 77, for N = 15, 150, 1500 (A,B, 
and C, respectively). The clusters applied where grown with 
a) f] = 2 and b) 77 = 4, consistent with the extreme val- 
ues of the maximum likelihood. Note that the products were 
normalized by their maximum value and that every second 
growth step was used only. The latter was done to avoid the 
products becoming too small. 



pronounced and the 77 value used in the simulations is eas- 
ily recovered. These products confirms that the number 
of Monte Carlo samples used in Eq. ([5]) are appropriate 
and more importantly that the method is directly appli- 
cable to experimental data for estimating an effective 77 
value or more generally the boundary condition function 
determining the growth rate. 



IV. DIMENSION AND a min 

The dimension of a cluster grown by this conformal map- 
ping technique is determined by the first term in the 
Laurant expansion of F[ n \ which will scale like 

~ n 1 / D \f\o @. The dimension is thus estimated 
by a direct fit of this scaling law as demonstrated in Fig. 
|3] for a cluster 80000 particles and 77 = 4.0. Using the 
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FIG. 3: (Color online) First Laurant coefficient versus cluster 
size for 77 = 4 and 80000 particles (with k = 2, see text). The 
added line is a fit of the fractal dimension D = 1.10 

conformal mapping technique we have grown clusters up 
to sizes 80000 particles with varying values of 77 in the 
interval 77 € [1,5]. Fig. [5] shows the results for the value 
of the dimension versus 77. As is clear from the figure, 
the value of the dimension decreases smoothly with 77, 
from the DLA value D = 1.71 for 77 = 1 down towards 
D ~ 1 for 77 — > 00. Hastings [8| presented arguments in 
favor of an upper critical dimension 7/ c = 4 for which the 
clusters become one-dimensional. We however do not ob- 
serve indications of this transition. As seen in Fig. 0] it 
is quite clear that the data smoothly bends away before 
reaching the point (77, D) = (4, 1) and only approaching 
the one-dimensional growth in the limit of large 77- values. 
We thus conclude that there do not exist a critical point 
at a finite 77. 

Halsey has computed a first-order correction to D 
for 77 < 4, obtaining D = 1 + |(4 — 77) + 0(4 - r/) 2 . 
This relation predicts a linear variation of slope \ around 
77 c = 4. As seen in Fig. 0]we do not observe this behavior. 
It is well know that the growth measure of a DBM 
model exhibits multifractal properties with a spectrum 
of growth exponents measured by local Hoelder expo- 
nents a [5| . The points of highest growth measures are 
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FIG. 4: (Color online) Computed values of the fractal dimen- 
sion versus r\. The dimension was extracted from a fit of the 
first Laurant coefficient Fi using clusters of sizes 20000-40000 
particles for r\ < 3.5 and 40000-80000 particles for r\ > 3.5. 
Each data point is averaged over 20 clusters and the error 
bars are estimated by the standard deviation. 



characterized by the minimum a- value, a m in- We have 
earlier determined this value using the iterated conformal 
mapping technique [ll| and extend it here to the DBM 
model. In this method, it is very easy to keep track 
of where the maximum growth probability is located as 
more particles are added. Let us assume that at the (n- 
l)'th growth step the site with the largest probability is 
located at the angle 9 max on the unit circle, i.e. for all 9 
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|$(n-l)'( e i0 m<M )| - |$(™-l)'(e i9 )| 



(8) 



When we add a new bump in the n'th growth step the 
position of maximal probability may not change (up to 
reparameterization of the angle 9 max ), or move to the 
new bump. We can easily find the reparameterized angle 
and determine the new position from 



max,n 

p 1 = max ■ 



1 



1 



l* (n) W.fl> Wm "))l 



|$W(e«»)| 



(9) 

If p™ ax ' n is located at 9 n we put 9 max — 9 n in the (n + 
l)'th growth step. Using conformal mappings, we have 
also previously estimated the critical branching angle as 
a function of 77 in the DBM model [l2j • 
Fig. [5]shows the results of a m i n vs. 77 and we observe that 
drain decreases from the DLA values a m i n = 0.68 down 
to a m i n = 0.5. It is obvious that a m i n = 0.5 corresponds 
to the Hocldcr exponent for a line. In consistency with 
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FIG. 5: (Color online) Computed values of a m i n versus r/ 
using the same clusters as in Fig. [4] 
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FIG. 6: (Color online) a m ; n versus the fractal dimension D. 
The data are the same as used in Figs. [4] and [5] 



the results in Fig. [4] we observe that the curve bends 
smoothly and that the one-dimensional growth is only 
obtained in the limit r\ — * 00. The last figure, Fig. [6l 
shows a m i n plotted vs. D. By extrapolation (as indi- 
cated by the line) we see that a m i n assumes its minimal 
value 0.5 at a dimension D = 1.0. 
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V. CONCLUSIONS 

The conclusions of our paper are twofold. Firstly, we 
have presented a method to extract the effective value 
of the growth exponent 77, for a time series of a growing 
aggregates, assuming an underlying mechanism based on 
the Dielectric Breakdown Model (DBM) model. The es- 
timate is based on a maximum likelihood method and 
converges rather well for the numerical data presented 
here. We believe this method should be directly appli- 
cable to experimental data when it is possible to extract 
intermediate steps in the formation of the aggregates. 
We urge the method to be used in for example viscous 
fingering experiments in random media 13|. Secondly, 
we have thoroughly investigated the scaling structure of 
DBM clusters as a function of the growth exponent t). 
Based on extensive numerical simulations wc do not find 



support for the conjecture that the growth becomes one- 
dimensional at the critical value r\ c = 4 [3, @]. On the 
contrary, our results indicate that there do not exist a 
critical point for at finite 77-value and that the scaling 
exponent of the maximal growth site a m i n assumes its 
minimal value 0.5 when the growth becomes non-fractal. 
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